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Abstract 

We report intimate relations between topological properties of full gapped spin-triplet super- 
conductors with time-reversal invariance and the Fermi surface topology in the normal states. An 
efficient method to calculate the Z2 invariants and the winding number for the spin-triplet su- 
perconductors is developed, and connections between these topological invariants and the Fermi 
surface structures in the normal states are pointed out. We also obtain a correspondence be- 
tween the Fermi surface topology and gapless surface states in the superconducting states. The 
correspondence is inherent to spin-triplet superconductivity. 

PACS numbers: 



Electronic address: msato@issp.u-tokyo.ac.jp 



Search for possible states of quantum matter is one of the central issues in condensed 
matter physics. In addition to local order parameters, gapped states can be characterized 
by topological invariants which are constructed from the wave functions. The quantum Hall 
state is a prominent example of such topological states, in which the Hall conductance is 
identified with the topological number introduced by Thouless, Kohmoto, Nightingale and 
den Nijs (TKNN) While the TKNN number is non-trivial only when the time-reversal 
symmetry is absent, recently new topological invariants, namely the Z2 invariants, were 
introduced in order to distinguish a topological state with time-reversal invariance from 
ordinary band insulators . For such a "topological insulator', the bulk-edge 

correspondence between the topological invariants in the bulk and gapless edge (or surface) 

n n Q 

states on the boundary was discussed in a similar manner to the quantum Hall state |2J, |5|, |9(. 
The topologically protected gapless state is an origin of dissipationless (spin) Hall effects, 



which inspire an application to spintronics lOl. 111!]. 



In this paper, using the topological invariants, we study topological properties of another 
class of gapped systems, spin-triplet superconductors. Although conventional s-wave super- 
conductors are topologically trivial, it is known that an unconventional superconductor can 
be topologically non-trivial 12J, [13j, ll4| • In the following, we will develop a powerful method 
to evaluate the topological invariants for spin-triplet superconductors, and find an intimate 
relation between the topological properties in the spin-triplet superconducting state and 
those in the normal state. In particular, from topological arguments based on the bulk-edge 
correspondence, we derive formulas between the gapless surface (edge) state on the bound- 
ary of the three dimensional (3D) (two dimensional (2D)) spin-triplet superconductor and 
the topological invariants of the Fermi surface in the normal state. Although the number 
iVo of the gapless surface (or edge) states itself depends on the details of the gap function, 
it will be shown that the index (— 1) does not depend on them and is directly related to 
the Fermi surface topology in the normal state. We also introduce a tight-binding lattice 
model, and confirm the results by numerical calculations. 

In the following, we consider mainly full gapped spin-triplet superconductors with time- 
reversal invariance. A generalization to those without time-reversal invariance will be men- 
tioned in the last part of this paper briefly. 

Let us start with the single-band description of a spin-triplet superconducting state with 
time-reversal invariance. (Generalization to the multi-band description is presented later.) 
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The Hamiltonian H of a spin-triplet superconductor in a single-band is given by 

n = E £ (*)4^ + \ E ( A -' ( fe ) c ^c V + h - c ) ' W 

fc,cr k,a,cr' 

where (cfc j<7 ) denotes a creation (annihilation) operator of the electron, e(fe) the dis- 
persion of the electron in the normal state and A(fe) the gap function given by A(fc) = 
id(k) ■ cr<7 2 . d(k) are odd functions and cr are the Pauli matrices. By rewritten H as 

H = I E 4# 4 = (4 iCT , c _ fc)(7 ), (2) 

it is found that the spin-triplet superconducting state is describe by the 4x4 Bogoliubov-de 
Gennes (BdG) Hamiltonian, 

Urn^ A(») \ 

\ A(fc)» -e(fc)l 2x2 j 

We assume that the normal state has the inversion symmetry and the time-reversal invari- 
ance, so e(— k) = e(k). From the time-reversal invariance of H(k) 

0iJ(fc)0- 1 = H(-k)*, = f ^ ° ) , (4) 

\ io-2 J 

d{k) should be real. Eigenstates of H{k) with negative energies E{k) < are occupied in 
the ground state of the superconducting state. 

The essential ingredient of our argument is the following "symmetry" for the spin-triplet 
superconductor. Since the parity of the gap function is odd, d(—k) = —d(k), the BdG 
Hamiltonian of the spin-triplet superconductor has the symmetry 

UH{k)W = H(-k), n 2 = 1 (5) 

with 

II I 12X2 | =1 2 X2®73. (6) 
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Using this, we will study topological properties for the spin triplet superconductor. 

Let us consider special points k = T a in the Brillouin zone which are time-reversal 
invariant and satisfy — r a = r a + G for a reciprocal-lattice vector G. In terms of the 



primitive reciprocal lattice vectors bj, the time-reversal invariant momenta T a are expressed 
as 



r a =(ni,n 2 ) = ( n ibi + n 2 b 2 )/2 for two dimensions, (7) 
r o =(ni,n2,n 3 ) = {n\bi + n 2 b 2 + n 3 b 3 )/2, for three dimensions, (8) 

with rij = 0, 1. At these momenta, the time-reversal invariance (jlj) reduces to 0-ff(r a )@ _1 = 
H(T a )*, since H{k) satisfies H{k + G) = H{k). This implies that an occupied eigenstate 
\uniXa)) {n — 1)2) has the same energy as its Kramers partner Q\u n (T a ))* . In addition, 
from the additional symmetry (151) , we have [if(T a ),n] = 0. So the Kramers doublet of the 
occupied states has the same eigenvalue of II. The eigenvalue of II is given by 

vr a = -sgne(r a ) (9) 

since H(T a ) = e(T a )Tl. 

The eigenvalues {vr a } have the following interesting properties: a) They are defined only 
at the time-reversal momenta {T a }. b) They only take 7r a = ±1. c) Their values can 
change only when the gap of the system closes. To see the last property c), consider the 
quasiparticle spectrum, E{k) = ±^/e(k) 2 + d(k) 2 , which is obtained by diagonalizing H(k). 
The gap of the system 2\E(k)\ closes when e(k) = d(k) = 0. At the time-reversal momenta, 
the d vector vanishes identically, d(T a ) = 0, so only e(T a ) = is required for gap closing. 
Therefore, the gap closes when n a changes. 

The above properties suggest a connection between the Z 2 invariants introduced in 
and {vr a }: The Z 2 numbers are calculated from the quantities {5 a } 

vgetRrg] 

a PfMr a )] ' ^ 

where w(T a ) nm is the anti-symmetric U(2) matrix connecting the occupied states \u n (T a )) 
{n = 1,2) with their Kramers partners 9|w n (r a ))*, w(T a ) nm = (w n (r a )|6|w m (r a ))*, and Pf 
denotes its Pfaflian. While the quantities {5 a } depend on the gauge (or phase choice) of 
the occupied states, their gauge-independent combinations define the Z 2 invariants, v for 
two dimensions and (/i = 1,2,3,0) for three dimensions (|: (— l) v = Yin =o l ^a=(ru,n 2 ) 
for two dimensions, and {-l) v ° = H n . =01 5 a={ni>n2>n3) , (-l) Uk = Un^ k =o,i;n k =i S a={n u n 2 ,n 3 ) 
(k = 1,2,3) for three dimensions. We notice here that the quantities {5 a } have properties 
similar to those of {vr a }: A) They are defined only at the time-reversal invariant momenta. 



B) They only take 5 a = ±1 since Pf[u>(T a )] 2 = det[w(r a )]. C) With fixing the gauge (or 
phase choice) of the occupied states, their values can change only when the gap of the system 
closes. The last property C) is obvious because their gauge independent combinations v and 
can change only when the gap of the system closes. 

These similarities suggest that the relation S a = ir a holds with a suitable phase choice 
of the occupied states. Indeed, we can prove it by using a similar technique developed in 
{(J with the replacement of the inversion symmetry P by the symmetry II in the argument. 
As a result, we obtain useful formulas of the Z 2 invariants for the time- reversal invariant 
spin-triplet superconductor, 

{-If = Yl s g ne ( r a=(ni,n 2 )) 3 for two dimensions, (11) 

n,=0,l 

(~ i y° = I [ s g ne ( r a=(ni,n 2 ,n 3 )), (-1)"* = JI sgne(r a=(ni ,„ 2 ,„ 3 )) , for three dimensions. 

«j=0,l nj^ k =0,l;n k =l 

(12) 

Note that the Z 2 numbers v and are mod 2 integers, which are identified with v + 2 and 
z/^ + 2, respectively. Thus the Z 2 numbers are non-trivial (trivial) when they are odd (even). 

Here we find that the right-hand sides of Eqs. fllip and (1121) have their own topological 
meanings related to the Fermi surface structure in the normal state: For Eq.( TlT]) . by using 
the relation e(— k) = e(fe), it is found that 

J] sgne(r a=(ni , n2) ) = (-1)»<*W (13) 

n,j=0,l 

where po(Sp) is the number of different connected components of the Fermi surface in the 
normal state. Also for Eq ffl2l) . we obtain 

J] sgne(r a=(ni , n2 , n3) ) = (-1)*<W, 

n J= 0,l 

J] sgn6(r a=(niin2 , n3) ) = (-l)^(^), (14) 

«j^ft=0,l;nfc=l 

where x{Sf) is the Euler characteristic of the Fermi surface, and Po(Cfc) is the number of 
different connected components of the intersection Ck between the Fermi surface and the 
time-reversal invariant plane with k = bk/2. (For a single connected Fermi surface, the 
Euler characteristic is given by x{Sf) — 2(1 — g) with g the genus of the Fermi surface. 
When there are multiple connected components of the Fermi surface, x{^f) is the sum of 



the Euler characteristics of each component.) We illustrate pq(Sf), x(Sf) and po(Cfc) hi 
Figs{T] and [2j Eqs. ffTBl and (fT4l) are confirmed by these examples. These quantities, po(Sf), 
x(Sp) and po(Cfc) are topological invariants of the Fermi surface, and they do not change the 
values under deformations of the Fermi surface unless the Fermi surface crosses one of the 
time-reversal invariant momenta. Therefore, Eqs. fllip - CH]) make connections between the 
topological invariants in two different phases, i. e, . the Z2 invariants in the superconducting 
phase and po(Sf), x(Sf) and po(Cfc) in the normal phase: 

(-If = (-1)po( s p), for two dimensions, (15) 
(_rp = (_i)x(Sf)/2 ) ^_ x y k = (_ 1 )po(o*) j for three dimensions, (16) 

An important physical consequence of our formulas ( fTBl) and ( fT6|) is that one can obtain 
useful information about gapless surface (or edge) states in the spin-triplet superconductor 
from the knowledge of the Fermi surface topology. From the bulk-edge correspondence, 
a non-trivial Z2 number of a bulk gapped system implies the existence of a gapless state 
localized on the boundary [5J. For time- reversal invariant systems in two dimensions, the 

n n 

gapless state is non-chiral and its Kramers doublet forms a helical pair [11|, [15J. The helical 
edge pair also satisfies the Majorana condition in the present case, because of the particle- 
hole symmetry of the superconducting system. From a topological argument similar to that 
in {5I, it is shown that an odd (even) number of gapless helical Majorana pairs exist on 
each edge when (— l) u = —1 ((— l) u = 1). Thus from (fT5]) . we find the following connection 
between the number N of the gapless helical Majorana pairs on each edge and the topological 
invariant po(Sp) of the Fermi surface, 

(-1)*° = (17) 

This formula implies that when po(Sp) is odd, iVo cannot be zero, and at least one gapless 
helical Majorana state should exist on each edge. 

For 3D time-reversal invariant spin-triplet superconductors, the gapless boundary state 
is a 2D massless Majorana fermion. By generalizing the argument in {5] to this case, we 
have the following two properties of the surface state. 1) The number No of 2D gapless 
Majorana fermions on a boundary surface is related to the topological number vq by the 
equation (— 1)^° = (— l) u °. 2) When (— l) u ° = 1, a non-trivial z/j implies the existence of 2D 
gapless Majorana fermions on surfaces determined by Uf To specify the surfaces, consider a 



surface G which is perpendicular to a reciprocal-lattice vector G. If the surface G satisfies 
G 7^ Ylii^i + 2mj)6i for any integers rrii, then there exist 2D gapless Majorana fermions 
on the surface. Combining the former property with (j!6p . we have a relation between the 
gapless surface state of a 3D time-reversal invariant spin-triplet superconductor and its Fermi 
surface topology as, 

(-1)^ = (_ 1 )x(Sf)/2 ) (lg) 

where N Q the number of the 2D gapless Majorana fermions on the boundary surface. More- 
over, taking into account the latter property as well, we obtain the following predictions, (i) 
When the Fermi surface satisfies (— l)^ 5 ^)/ 2 = —1, an odd number of 2D gapless Majorana 
fermions exist on each boundary surface. In particular, at least one gapless Majorana fermion 
exists on each boundary surface, (ii) When the Fermi surface satisfies (— l)^ 5 ^)/ 2 = 1, 
the number of the 2D gapless Majorana fermions on a boundary surface is even. Then if 
the surface G satisfies G ^ J2i(Po(Ci) + 2m^)bi with arbitrary integers at least two 
2D massless Majorana fermions exist on the boundary surface G. On the other hand, if 
G = ^2i{po{Ci) + 2mi)bi with integers rrii, no gapless Majorana fermion is possible on the 
surface G. 

In Table (TJ we summarize the relations between the Fermi surface topology and the 

n 

boundary gapless state [24j . Later, we will check these results by using concrete models. 



For 3D time-reversal invariant spin-triplet superconductors, it is also known that there 
exists another topological invariant z/ w called the winding number [if], Now we will 
derive a useful formula for z/ w and show that z/ w also has an intimate relation to the Fermi 
surface topology. In the single band description, the winding number i/ w is given by 



UW '~ 127T 2 



-L / dkW^Tjed&djfjcdkVd, (19) 

2ir z J T 3 

where T 3 denotes the first Brillouin zone, and fj a (k) = r] a (k) / ^Jr] a {k) 2 with rj a (k) = 
(d(k) , e(fc)) . v w counts the number of times the unit vector fj a wraps the 3D sphere S 3 
(f)a = 1) when we sweep T 3 . In order for fj a to wind S 3 , it is necessary to pass the poles of 
S 3 defined by rj = {i]i,i]2,i] 3 ) = d = 0. So consider the set of zeros k* satisfying r)(k*) = 0. 
From the topological nature of z/ w , we can rescale e(fc) as e(fc) — > ae(k) (a 1) without 
changing the value of z/ w . Then it is found that only neighborhoods of the zeros contribute 
to z/ w if a is small enough. By expanding r\ a as r\i = djdi(k*)(k — k*)j + ■ • ■ , (z = 1, 2, 3), 



a) 2D case 



(_1)Po(Sf) = _1 


JV = 1,3,5,--- 


(_l)Po(S F ) = x 


N = 0,2,4,- •• 



b) 3D case 



(_1)x(Sf)/2 = _1 


Ao = 1,3,5, ••• 


(_1)x(Sf)/2 = 1 


On a surface G = ^2i(po(Ci) + 2rrn)b{ 


iV = 0,2,4,- • • 


On a surface G ^ Sj(Po(Q) + 2m,j)&j 


iV = 2,4,6,- •• 



TABLE I: Topological invariants of the Fermi surface and the possible number Ao of gapless 
boundary states for full gapped time-reversal invariant spin-triplet superconductors, a) 2D case. 
Here po(Sp) denotes the number of connected components of the Fermi surface, and Ao the possible 
number of gapless helical Majorana pairs on an edge, b) 3D case. Here x(S'f) is the Euler 
characteristic of the Fermi surface, po(Cfc) the number of different connected components of the 
intersection C\~ between the Fermi surface and the time-reversal invariant plane with k = &&/2, 
and nii integers. The surface G is perpendicular to the reciprocal-lattice vector G, and Ao is the 
possible number of 2D Majorana fermion on the surface G. 

7/4 = e(k*) (<C 1), the contribution from the zero k* is evaluated as 

u w (k*) = -isgn(e(fc*))sgn (det^-t^fe*))) . (20) 

(When det (djdi(k*)) = 0, Eq.(|20l is generalized to 

u w (k*) = -hgn(e(k))i(k*), (21) 

where i(k*) denotes the Poincare-Hopf index [l8j of the zero k*.) Summing up the contri- 
butions of all zeros, we have 

?7(fc*)=0 

From ( 122|) . we can show that z/ w is also related to x(Sf)- For simplicity, suppose that the set 
of zeros k* contains only the time- reversal invariant points {T a }. (T a is always a zero since it 
satisfies d(T a ) = 0.) Dividing the set of zeros into two subsets, T± = {r a ; sgne(r a ) = ±1}, 
we obtain v w = — J2r a er+ i(T a )/2 + X]r a er_ *(r a )/2. Then by using the Poincare-Hopf 
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theorem ^2 k * i(k*) = [19J], it is recast into z/ w = X]r a er_ i(Fa)- Here i(T a ) is an odd 
integer because of d(—k) = —d(k). Therefore, z/ w is an odd (even) integer if T_ has an odd 
(even) number of elements. From this, we obtain the relation 

(_l)-w = Yl sgne(r a=(niin2 , n3) ). (23) 

nj=0,l 

Combining this with (THl) and (fl8l) . we find that z/ w is also related to the Euler characteristic 
x(Sf) and the number N of 2D gapless surface states as 

(_l)-» = (_!)x(S P )/2 = (_i)JVb_ (24) 

Let us now illustrate our results with simple and important examples. In FigJTJ we 
illustrate possible Fermi surfaces in the normal state and the corresponding po(Sp) in two 
dimensions. We also present the energy spectra for the corresponding superconducting states 
with edges. To obtain the energy spectra, we use the lattice model of the superconducting 
state with d = cf(sin k x x + sin k y y), 

1 * — \ i / ta Ida ■ <ra<2 \ 

W = ;E C ^ H V = 3 ), (25) 

z ij \ -i<Ta 2 dji -Uj J 

where c] = (ct CT ,c i(J ), and Uj and d y are given by % = -t x (5^ j+A + 8 jM& ) - t y (5 i)j+y + 
$j,i+y) ~ ^iji {dx)ij = ~ i(d/2)(5jj+x ~ ^i,j+x), (d y )ij = —i( K d/2)(5j^ + y — 5ij + y), (d z )ij = 0. 
The spectra are calculated for the system with two edges at % x = 0, 50 under the periodic 
boundary condition in the y-direction. In FigJTJ k y denotes the momentum in the ^-direction. 
While no gapless edge state exists in FigfjTJa), it is found that there exist gapless edge states 
in the bulk gap in Figs. [T]b) and c). The relation Eq. fll7l) holds in FigJTJ 

In FigfJl we show various Fermi surfaces in the first Brillouin zone and their topological 
numbers, x(Sf) andpo(C?) {i = 1,2,3) in three dimensions. In addition, we present gapless 
2D Majorana surface states for the the superconducting states with d{k) = sin k x x+sin k y y+ 
sin k z z. This figure also confirms the connection between the gapless surface states and 
the Fermi surface topology: The relation (— 1)x( s f)/ 2 = (— 1)^° holds for all the cases. 
Furthermore, in the cases with (— l)*^)/ 2 = 1 [i.e.. FigfS] b) and d) ), there exist a non- 
zero even number of 2D gapless Majorana fermions on a surface G ^ Y2i(Po(^i) + 2mj6j) 
with integers m^. Note that in Figj2j only the 001 surface in Figj2]b) does not satisfy this 
condition. In this case, we have G = b 3 , and it coincides with J2iPo{Ci)bi = b 3 . From (1221) . 



a) k y 
n 
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FIG. 1: (color online). The Fermi surfaces in the normal state and the edge states in 2D time- 
reversal invariant spin-triplet superconducting state. (Top row) The Fermi surfaces and 7r a = ± 
at the time-reversal invariant momenta. po(= Pq(Sf)) is the number of the connected components 
of the Fermi surface. (Bottom row) The energy spectra of the corresponding superconducting 
states described by ([25]) with edges at i x = and i x = 50. Here k y denotes the momentum in the 
y-direction, and Nq the number of gapless helical edge states. We set the parameters of the lattice 
model (f25|) as (a) t x = 0.4, t y = 1, \x = —1, d = 0.5, (b) t x = t y = 1, /i = —1, <i = 0.5, and (c) 
t x = 1, t y = 0.4, fi = —1, d = 0.5, respectively. 



we find that u w 's for this gap function are (a) z/ w = 1, (b) v w = 0, (c) z/ w = —1, and (d) 
u w = —2, respectively. These values are also consistent with Eq. (l24l) . 

So far we have considered the single-band superconductor. However, the formulas pip 
and (fl~2]) can be generalized to multi-band systems. To see this, consider a multi-band 
system which has the inversion symmetry and the time-reversal invariance in the normal 
state. If we assume that the parity operator transforms only the momentum as k — ► —k 



251 ]. then the Hamiltonian in the normal state is given by a 2N x 2N matrix S(k) satisfying 



S(—k) = £{k). (N is the number of the bands.) Odd-parity superconducting states for this 
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FIG. 2: (color online). Various Fermi surface topologies in three dimensions and the corresponding 
gapless surface states for the 3D time-reversal invariant spin-triplet superconductor with d(k) = 
sin k x x + sin k y y + sin k z z. (Top row) Fermi surfaces in the first Brillouin zone and their topological 
invariants, x(Sf)',\po(Ci),po(C2),Po(C3)]. The green circles are Cj (j = 1,2,3) for d). (Middle 
and bottom rows) The corresponding surface states on the Brillouin zone for 001 surface (middle) 
and 100 surface (bottom) in the superconducting state. The blue solid circles symbolize the Dirac 
cones of the 2D gapless Majorana fermions, and the energies of the surface states become zero at 
the time-reversal invariant momenta enclosed by the blue circles. Nq denotes the number of 2D 
gapless Majorana states on each surface. 



system are described by the generalized BdG Hamiltonian 

m) - ( m A(k) ) , («, 
\A( k y -£(k)J 

where the gap function A(k) is a 2N x 2N matrix with odd parity, A(— k) = —A(k). H(k) 
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has the property 

n//(fc)n t = H(-k), n 2 = 1 (27) 

with II = l2AT X 2Af <8> t 3 , and for A; = T a , -H"(fe) becomes -H"(r a ) = S(T a ) ® r 3 , thus in a similar 
manner as the single-band case, it is shown that 

N 

(m,n 2 )))> f° r two dimensions, (28) 

nj=0,l m=l 
iV 

(~ 1)"° = n n s g n (^2m(r a= ( ni , n2in3 ))), 

nj=0,l m=l 

TV 

(-1)"* If [J sgn(E 2m (r a=(niin2in3 ))), for three dimensions, (29) 

nj-t k =Q,l;n k =l m=l 

where E n (T a ) (n = 1,---2N) are the eigenvalues of £{k) at fe = r a , and we have set 
-E^m^a) = -E l 2m^i(r a ) using the Kramers degeneracy. For a filled or empty band in the 
normal state, the signatures of E n (T a ) are the same for all the time- reversal points, so their 
contributions to fl28|) and fl29|) are canceled. Therefore, in order to evaluate the Zi numbers, 
it is enough to consider bands with the Fermi surfaces. Again it is evident that topological 
properties of the spin-triplet superconducting state are closely related to the topology of the 
Fermi surface. 

Finally we make several comments in order, a) Although we have assumed that the 
normal state has the inversion symmetry, our formulas (1TT1) and (1T21) (or (1251) and (1291) ) 
could be useful even for the systems which do not have the inversion symmetry in the normal 
state: Adiabatic continuity allows us to calculate the topological invariants if the system 
is adiabatically connected to materials which have the inversion symmetry in the normal 
state. The topological invariants for a class of noncentrosymmetric superconductors can be 



calculated in this manner 



20 



2l|. b) For spin-singlet superconductors, due to the inversion 



symmetry, their Z 2 numbers are calculated by the technique developed in 



is found that all the Zo numbers are trivial 



However, it 



22]. Therefore, the correspondence between the 



Fermi surface topology and the gapless surface state discussed in this paper are inherent to 
spin-triplet superconductors, c) In this paper we have focused on the time-reversal invariant 
spin-triplet superconductors. Here we mention a generalization to the time-reversal breaking 
case in brief. For 2D chiral spin-triplet superconductors such as a p+ip state, the topological 
properties are determined by the TKNN number ^tknn- I n a similar manner to u w , in the 

12 



single-band description, it can be shown that the TKNN number is related to the Fermi 
surface topology by the equation 



(30) 



where Po(Sf) is the number of the connected components of the Fermi surface 22]. This 



relation gives a simple explanation of the quantum phase transition from the weak paring 



phase to the strong one discussed in 23]. This phase transition is accompanied with dis- 
appearance of the Fermi surface, thus Po(Sf) — 1 — > Po{Sf) = 0. From the above relation, 
this cause a change of z/ TK nn, which brings about different topological properties between 
the weak and strong phases. 

The author would like to thank Y. Asano, Y. Tanaka, X.-L. Qi, R. Roy, and Y.-S. Wu 
for helpful discussions. 
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